Incomplete current fluctuation theorems for a four-terminal model 
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We demonstrate the validity of the current fluctuation theorem for a double quantum dot sur- 
rounded by four terminals within the Born-, Markov- and secular approximations beyond the 
Coulomb-blockade regime. The electronic tunneling to two fermionic contacts conserves the total 
number of electrons, and the internal tunneling is phonon-assisted by two bosonic baths. Adapted 
choice of thermodynamic parameters between the baths may drive a current against an existing 
electric or thermal gradient. We study the apparent violation of the fluctuation theorem when only 
some of the energy and matter currents are monitored. 
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Fluctuation theorems (FTs) connect forward and back- 
ward probabilities for processes associated with a definite 
exchange of entropy 1 . Thereby they relate rather sophis- 
ticated and hard-to-calculate quantities with simple and 
universal thermodynamic parameters, which constitutes 
part of their attractiveness. When the matter and en- 
ergy currents are tracked over a certain period of time, 
there exist simple versions of the current (or Full Count- 
ing Statistics) FT 2-5 . Significant progress made in the 
monitoring of electronic tunneling events through quan- 
tum dots (QDs) 6 ' 7 by using capacitively coupled quan- 
tum point contacts (QPCs) has led to an accurate under- 
standing of Full Counting Statistics 8 . Unfortunately, the 
QPC signal originating from a monitored single dot does 
not allow to reconstruct bi-directional tunneling events. 
Therefore, monitored double quantum dots allowing for 
bi-directional counting have entered the focus of inter- 
est 9,10 and appear as ideal testbeds to check current FTs. 



It is therefore essential to identify processes that may 
lead to modifications of the FT in an experimental setup. 
The FT has been found to be modified due to true quan- 
tum effects such as Berry phases 11 or interference ef- 
fccts 12 - 14 . It may however also be modified due to de- 
tector back-action - e.g., the back-action of QPC on a 
monitored double quantum dot 15,16 or the influence of a 
single electron transistor monitoring another one 17,18 - 
or simply when detailed balance is explicitly broken via 
feedback control 19 . 



Here, we will argue at the example of phonon-assisted 
tunneling that an apparent violation of the current FT 
may also arise due to ignored couplings with further 
baths (that might resist an experimental monitoring). 



I. MODEL AND METHODS 
A. Hamiltonian 

We consider a double quantum dot system (see Fig. 1) 

Hs = CAd^dA + Z-Bd^dB + Ud A dxd^dB (1) 

with on-site energies £a and £b and Coulomb-interaction 
U. Without loss of generality (mirror symmetry) we as- 
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FIG. 1: (Color online) Two quantum dots A and B are tunnel- 
coupled to adjacent fermionic contacts described by Fermi- 
distributions }l and fa with rates Fl and Fr, respectively. 
Electronic transitions between the QDs are only possible in 
the singly-charged sector via spontaneous emission or absorp- 
tion of a boson with energy transfer £b — £a from upper or 
lower phonon baths (defined by Bose-Einstein distributions 
nu and ud) with rates Y\j and Td. We explicitly allow for 
double occupation of the system, such that the coupling be- 
tween particles and energy transferred to the fermionic reser- 
voirs is not tight: Depending on whether the inert QD is 
occupied or not, the energy exchanged (red, in brackets) with 
the fermionic baths may vary. The use of multiple particle 
counting fields specific to both reservoir and transferred en- 
ergy facilitates the calculation of the complete particle-energy 
counting statistics. 



This paper is organized as follows: In Sec. I we in- 
troduce our model and the method, followed by a veri- 
fication of the multi-terminal FT in Sec. II. We discuss 
how the multi-terminal FT is modified when only partial 
information is available in Sec. III. 



sume that £a < £b- The system is surrounded by two 
bosonic a S {U, D} and two fermionic a £ {L, R} baths 

Kb = ^ £fc, Q c^ Q c fciQ + '^2^k,abl a b k ,a , (2) 

k , at k,a 
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that arc assumed to remain in thermal equilibrium 
throughout. Each QD is coupled to its adjacent fermionic 
contact by the tunneling Hamiltonian 



Hf = d A ^2 ^,l4,l + rf B ^2 tk,Rc{ R + h.. 



(3) 



where the tunneling amplitudes tk, a lead to effective tun- 
neling rates r a = T a (uj) = 2-7r^ fc \tk, a I ~ £fc,a) that 
we assume to be energy independent (wideband limit). 
In contrast, the transition A o B is phonon-assisted via 
the interaction 



nf = (d A 4 
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such that - under the rotating wave approximation - an 
electron jump between the QDs goes with either emission 
or absorption of a boson from upper or lower bath 20 . We 
again summarize the corresponding amplitudes in effec- 
tive energy-independent phonon-assisted electron tunnel- 
ing rates T a = F CT (w) = 2nJ2k l^k^H^ — Uk,o)- The 
resulting total model is described by the sum of all Hamil- 
tonians H = Hs + ~Hb + Hf + H\ h and is visualized in 
Fig. 1. 



B. Liouvillian 

We assume to be in the sequential tunneling regime, 
where second order perturbation theory in the couplings 
tka and hka to the contacts is a good approximation. 
This regime can e.g. be achieved when all tunneling rates 
are small in comparison to the reservoir temperatures 
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More generally, Kondo physics is 



expected to be negligible when the reservoir tempera- 
tures are larger than the Kondo temperature 23 . In this 
regime, performing the Born, Markov, and secular ap- 
proximations 24 yields a master equation that is expected 
to yield valid results. It involves only the system density 
matrix, is of Lindblad-form 25 , and in the (localized) sys- 
tem energy eigenbasis (Hs |00) = 0, Hs [ 10) = £a 1 10), 
Us |01) = £ B |01), and H s |11) = (e A + £b + U) |11» it 
assumes the form of a simple rate equation as long as 
the energy levels of Hs are non-degenerate (recall that 
£b > £a)- In principle, the rates may be calculated rig- 
orously but the result is also evident from Fermis Golden 
rule. We assume the reservoirs to be in thermal equi- 
librium throughout, i.e., expectation values are given by 
the Fcrmi-Dirac 



(4,a C fc,«) = /«( e fc,a) 
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or Bosc-Einstein 
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(5) 



(6) 



distributions, respectively, where (3 a and /3 ff denote the 
inverse electronic or bosonic bath temperatures and /z Q 



and \i a the respective chemical potentials. Writing the 
diagonal entries of the system density matrix in a vector 
Ps = (poo,oo, pio,io, poifli, Pn,u), the master equation to 
this order assumes the form ps = Cps, where the Liou- 
villian superoperators are given by 
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Here, we have used the abbreviations 

n a = n a (e B - £a) , 

fh = /l(£a) i /r = /r(eb) , 

h = h(tA + U), fR = f R (e B + U), (8) 

compare Eqns. (5) and (6). In the Coulomb-blockade 
limit, the doubly occupied state always decays, which is 
formally expressed by the limits fa — > and fa — > 0, 
and the top-left 3x3 submatrix of the above equation 
reproduces previous models in the literature 26 . 

It should be noted that the above Liouvillian (7) only 
satisfies local detailed balance 27 : When the system is 
only coupled to a single junction S G {L, R,U, D}, 
detailed balance L^p\ = CJ^p\ - where p\ denotes 
the corresponding single-junction stationary state - is 
obeyed. When it is coupled to all four terminals how- 
ever, global detailed balance is broken. 



C. Conditional Master Equation 

Being responsible for transitions between different sys- 
tem states, the off-diagonal matrix elements of the Liou- 
villians in Eqns. (7) can be used to set up a connected 
set of equations for density matrices p^ n \t) conditioned 
upon the number of particles n that have tunneled to a 
respective reservoir 28 . That is, although the state of the 
system density matrix does neither yield information on 
the number of particles nor the energy transferred to a 
certain bath, the full time-record of a single trajectory 
would reveal this information. When we only focus on a 
single reservoir to and out of which we count all particle 
jumps, the conditioned density matrices follow an equa- 

+ £_ (0 ("+ 1 >, 



tion of the form = C p {n) + £+p (?1 
where £ + (£_) describes particle jumps into (out of) the 
reservoir of interest and Cq contains the remaining terms 
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(jumps to and from other reservoirs as well as the di- 
agonal matrix elements). Alternatively, such n-resolved 
master equations may also be obtained by performing 
the derivation after a detector has been added to the 
system 29 . By performing a Fourier transform 



ps(x,*) = £/4 n) (*) 



(9) 



one can convert the infinitely large n-resolved, condi- 
tional master equation to a four-dimensional one at the 
price of introducing the counting field \- Being e.g. inter- 
ested in the number of particles in the left contact, this 
formally corresponds to the replacements fi — > Jl& ix 
and (1 — /l) — > (1 — fL)e +%x in the off-diagonal matrix 
elements of Eqns. (7), respectively. 

However, here we are interested in the full energy- 
particle counting statistics, such that one has to treat e.g. 
single-electron jumps corresponding to an energy change 
of £a and to an energy change of £a + U in the system 
differently. Therefore, we introduce two counting fields 
for each fermionic contact and one for each bosonic con- 
tact. Note that in the Coulomb-blockade limit, energy 
and particle fluxes are tightly coupled (often also called 
strongly coupled 26 ), such that also for the fermionic con- 
tacts a single counting field would suffice. The Liouvillian 
as a function of energy-resolved counting fields reads 

= £(Xli,Xl2,Xri,XR2,Xui,Xdi) 
= £l (XLl , XL2 ) + Cr (XR1 1 XR2 ) 

+Mxui) + £ d (xdi), (io) 

where the explicit counting field dependence can be ob- 
tained from Eqns. (7) by performing the replacements 

/l -> / L e" iXL1 , (1 ~ h) -> (1 - A> +iXL1 , 
h -> / L e" iXL2 , (1 - h) -> (1 - / L )e +iXL2 , 



/r -> / R e- 1XR1 
h -> / R e- iXM 
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(1 - / R ) -> (1 - / R )e+ i ^ 1 , 

(1 - / R ) -> (1 - / R )e +i *** , 

(1 + nu) -> (l + nu)e +ixul , 

(l + n D ) -> (l + « D )e +i * D1 (ll) 

in the off-diagonal matrix elements of the Liouvillian. 
The counting fields are thereby related to specific transi- 
tion frequencies in the system 

XLl <">• £A , XL2 •!-> £a + U , 

Xri £b , XR2 O £ B + U , 

Xui ^ £b - £a , 

Xdi -> £b - £a ■ (12) 

Therefore, we can now check the FT for fermions and 
bosons not only separately 30 , but in a single model where 
both species interact. 

D. Full Counting Statistics 

From the Liouvillian (10), one obtains the full statistics 
in a specific reservoir-energy channel by taking deriva- 



tives with respect to the counting field, e.g. for the mo- 
ments 

(nUt)) = (~id Xl ) k M( Xl t)\ x=0 , (13) 
with the Moment-Generating-Function (MGF) 

M( X ,i)=^{e C{X)t p} , (14) 

where p is the stationary density matrix defined by 
£(0)p = 0. 

Most simple, the current as the time derivative of the 
first moment can be evaluated as (see e.g. 31 ) 



(i) = -m{d x L(x)\ x =op} . 



(15) 



It is generally much more difficult to reconstruct the full 
probability distribution from the MGF as this involves 
the inverse Fourier transform of the MGF - compare 
Eq. (9) 



7T 



(16) 



Provided a situation where Ai (x) is the only eigenvalue of 
the Liouvillian with Ai(0) = (but see e.g. 32 for a more 
generalized treatment) the cumulant-gcncrating function 
(CGF) becomes linear in time 7 



C( X ,t)=ln[M(x,t)} ->Ai(x)t, 



(17) 



in the long-term limit (which will be denoted by an arrow 
further-on). Therefore, X\(x) can be interpreted as the 
long-term CGF for the current. 



II. COMPLETE FLUCTUATION THEOREM 

The characteristic polynomial of the Liouvillian (10) 
is given by V(X,x) = dct[£(%) — Al] (not shown for 
brevity) . We have found that it fulfills the analytic prop- 
erty 



V(X,- X )=V(X,x + iA), 
where for all A the shift is given by 



/ ^l(ea-Ml) \ 
p h (t A + \J-p L ) 

/3r(£b - Mr) 
/3r(£b + U - pr) 
/?u(£b - £a - Pu) 
\ /3 d (£b - £a - Pb) J 



(18) 



(19) 



The characteristic polynomial can be decomposed as 

4 

"D{X,x) = I! ~~ ^i(x)] f° r au \ where Xt(x) denote 

»=i 

the four eigenvalues of Liouvillian (10). It therefore fol- 
lows that all eigenvalues must obey the same symmetry. 
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In particular, we have also for the dominant eigenvalue 
the CGF for the current - the symmetry relation 



Ai(-%) = Ai(x + iA). 



(20) 



with the same shift as in Eq. (19). 

The current FT is given by the ratio of forward and 
backward probabilities for particle and energy exchange 
with multiple baths in the long-term limit. For the 4- 
tcrminal model considered here (see Fig. 1), it follows 
from symmetry (20) by basic properties of the inverse 
Fourier transform (see, e.g. the appendixes of Ref. 33 for 
a more detailed discussion) and reads 



P+n(t) 
p_ n (t) 



lira 



e nLl/3 L (£A-A I L) e nL2^L(£A+U-/j L ) x 



X g™Rl/3R(£B-A'H)gnR2/3R(i;B+U-^R) 



xc 



xe 



"ui^u(£b-£a-A 1 u) 
iDifc(£s- £a— fin) 



(21) 



Using that transferred energies and particle numbers are 
related by 



"LI 



nu = nui 



riL2 , n R - 
no = "-di 



™R1 + «R2 , 



AE L = n L i£A + n L2 (e A + U) , 

AE R = n R i£ B + n R2 (e B + U) , (22) 

AEu = nui(£B-£A), AS D = 7i D i(£ B - £a) , 

we find that the result (21) is completely consistent with 
predictions in the literature 34 as one would expect for an 
effective rate equation satisfying local detailed balance. 



once one is only interested in e.g. the total number of 
bosons, the impact of the two reservoirs adds up to a 
hypothetical single reservoir at some average occupation, 
similar to previous findings 35 . This leads to an additional 
analytic property of the characteristic polynomial, which 
also transfers to the dominant eigenvalue 



Ai(xLi,XL2,XRi,XR2,X P h + i/3u(£B - £a - Mu), 
Xph + i^d(£b - £a - Md)) 

Ai(xli,XL2,Xri>XR2, 



Xph + i In 



1 + n 



,X P h + iln 



1 + n 



where the average bosonic occupation is simply given by 
the weighted sum 



ru«u(£B - £a) +r D ra D (£ B - £a) 
n = - — — , (25) 



ru + r D 



which allows one to define an average boson tem- 
perature at vanishing boson chemical potentials via 
A>h(£B - £a) = In 



Electronic Transfer FT 



III. INCOMPLETE FLUCTUATION 
THEOREMS 

It is evident that the numbers of particles counted at 
all junctions are not independent. The total number of 
electrons is conserved for example. This results in fur- 
ther analytic properties of the characteristic polynomial 
P(A,x), which transfer to the long-term CGF (17). We 
note here the identity 



Ai (xli , XL2 , XRl 7 XR2 j Xui, Xdi) 



Ai(xli 

XL2 "+ 

XRi 4 

XR2 4 

Xui 4 
Xdi 4 



si - s L , 



S2 - 

51 ■ 

52 ' 



«L, 
SR, 

sr, 
■ sr, 
- sr) 
(23) 



for arbitrary shifts s\, S2, sl, and sr. 

It is also obvious that the structure of the two bosonic 
Liouvillians is identical and that their dependence on the 
respective bath occupation is linear. This implies that 



We are interested in the joint probability that n c \ elec- 
trons have left the system at the right junction (i.e., 
Uri + ?ir,2 = 4-n e i) and n e \ electrons have entered the 
system at the right junction (i.e., Tin + "L2 = — ^ei)- 
This probability can in the long-time limit be evaluated 



PnAt) = 



E 



E 



E w) 



Tl R.l> rl R2 Tl Ll' Tl L2 L n Ul,«Dl 

nRi+nR2="ci "Li+nL2=— n e \ 

4^y ax > 

— 7T 

(26) 



where we have eliminated the integrals by using that 



En 1 



2n5(x)- In order to relate the backward prob- 



ability to the above equation, we consider in the limit 
where the electron temperatures are equal Pi = /3r = f3 c \ 
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the identities 



Ai(-Xi) -XL, ~Xr, ~Xr, 0,0) 
Ai (xl + i/?ei(£A - Ml),Xl + i/3 e i(£A + f - Ml), 
Xr + i/3d(£B - Mr), Xr + i/3ei(£B + ?7 - Mr), 
iAj( £ b - £a - mu),i/3d(£b - £a - Md)) 

Ai (xl + iAi(£A - Ml), Xl + i/3 e i(£A + E/" - Ml), 
Xfl + i/? e i(£B - Mr), Xr + i/?d(£B + U - Mr), 
i/3 P h(£B - £a), i/3 p h(£B - £a 

Ai (xl + i/3ei(£A - Ml + C^/2) + ^/3 p h(£B - £a), 

Xl + i/3 e i(£A - Ml + U/2) + ^/3 P h(£B - £a), 

XR + i/?ci(£B - MR + U/2) - ^/3 p h(£B - £a), 

Xfi + i/3ei(£B - Mr + U/2) - -/3 p h(£B - £a), 

0,0 



Eq. (20) 



Eq. (24) 



(27) 



Performing a shift of the integration variables in the de- 
nominator, we obtain the electronic transfer FT 



R 



+ ™cl 



(t) 



,(*) 



el [^el(eB-£A+V)-|8ph(£B-£A)] 



(28) 



where V = ^ — /iR denotes the conventional bias volt- 
age. That is, with incomplete information, the full fluc- 
tuation theorem acquires the form of a modified FT with 
a shift term in the exponential. Note also that the elec- 
tronic transfer fluctuation theorem does not depend on 
the Coulomb interaction term U. 

This has to be contrasted with modifications of the 
FT in the literature where one only observes a modified 
temperature 10,15,16 . Instead, the apparent violation of 
the FT effectively mimics one found for a Maxwell demon 
model 19 . 

The same FT as in Eq. (28) would be ob- 
tained if e.g., only the number of electrons leav- 
ing at the right junction was counted. This is 
a consequence of electron number conservation, 
which is formally expressed by the symmetry 
Ai(xi, Xl, Xr, Xr> 0, 0) = Ai(0,0,xn - Xl,Xr - Xl,0,0) 
- see also Eq. (23). 

The shift term in the FT has the interesting conse- 
quence that e.g. at zero bias voltage V = 0, a current 
may still be generated from left to right ((3 C \ > /3 p h) or 
vice versa (/3 e i < /? P h)- The conventional FT is repro- 
duced when /3 e \ = /3 p h. Such a transport-without-bias 
behavior may also be generated by introducing asymmet- 
ric energy-dependent tunneling rates 36 or when a bosonic 
bath couples directly to the QD occupation 37 . Also at fi- 
nite bias voltages, the device may perform work by trans- 
porting electrons against an existing potential gradient, 



Eq. (23) «2 




Pa 


< /3 P h 


- /3 pil AE = 


= 2 




= 200 


-/3 ph AE-- 


= 0.02 



-2 2 

bias voltage V/U 



FIG. 2: (Color Online) Renormalized electron current Ir — 
Iri + Ir2 from Eq. (15) versus bias voltage for different 
phonon bath temperatures. When the boson temperature 
exceeds the electron temperature, the current at zero bias 
voltage is positive (dashed red curve), whereas the opposite 
is true for larger boson than electron temperatures (dotted 
blue curve) - compare the zoomed inset. Parameters have 
been chosen as (3 c iAE = /3 c i(£b — £a) = 2, /3 e iU = 20, 

Mi = +v/2 = —hr, nu — fj-D — o, r L = r R = r v = r D . 



see Fig. 2. In our case, the required energy is provided 
by the temperature difference between the boson and 
fermion reservoirs. 



B. Bosonic transfer FT 



Evaluating the direct boson current between upper and 
lower baths requires to evaluate the full fluctuation the- 
orem at titji = +n p h and riDi = — %>h ~ disregarding 
the number of tunneled electrons at the other junctions. 
Formally, this corresponds to 



+ 7T 

p _> _L / P >i(0 J 0,0,0,xtr,xr.)t-inph(Xcr-Xu)rf2 v 

^n ph -> 47r2 / e a x ■ 



(29) 
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We now consider the identities (again for similar elec- 
tronic temperatures /3l = /3r = Pel only) 



Ai(0,0 > 0,0,-xt/,-Xr>) 
Ai (i/3 e i(£A - Ml), i/3 e i( £ A + U - Ml), 

i/3 c i(£ B - MR.))iAsi(£B + ^ _ Mr), 
+ i/3u( £ B - £a - Mu), 



Eq. (20) 



Xd + i/3d(£b - £a - Md 

Ai (0,0,0,0, 

Xu + 1Aj(£b - £a - Mu) - i/3 e i(£B - £a + V"), 
Xd + i/3d(£b - £A - Md) - i^d(£B - £a + V) 



Eq. (23) 



(30) 



Finally, this implies that the bosonic transfer FT (see 
also e.g. 38,39 ) 



g1pli[^u(£B-£A-fln)-&(£B-£A-f"D)] (31) 



is not affected by the electronic transport at all, which 
is due to the fact that for the boson-system coupling in 
our model, the coupling between transferred particles and 
energy is tight. 



Finally, this implies that the combined bosonic FT 

P+n pb (t) 



P-n ph (t) 



D «ph[/3ph(E-B-£A)-/3 c l(£B-£A + V)] ^g^-j 



is now affected by the electronic transport. We have also 
calculated the total bosonic emission rate and find qual- 
itative agreement with the FT, see Fig. 3. 



= 0,8 



0,4 





< /9 P h 







p el AE = 1.6 
/3 eA AE = 1.31 
/3 e iAE = 0.8 



-0,4 1- 

_i_ 



-2 2 

bias voltage V / U 



C. Combined bosonic FT 

When we do not differentiate in which of the bosonic 
baths phonons are counted n p h = nyi + tt-di, the proba- 
bility to count a given number of phonons is given by 

p v J_ f p Ai (0,0,0,0, x,x)t-i« P hX r 7 v 

— 7T 

(32) 

We now consider the identities (again for similar elec- 
tronic temperatures /?l = Ai = /3 c i only) 

Ai(0, 0,0,0, -*,-*) Eq '^ 0) 

Ai(i/3 c l(£A - yUL),i^ol(£A + U - /i L ), 

i/9 e i(£B - Mr), i/3 c i(£B + U - Mr), 

X + i/?u(£B - £a - Mu), 

■ n , \\ E q- (23) 

X + iPd(£b - £a - Md) J = 

Ai (0,0, 0,0, 

X + i/?u(£B - £A - Mu) - i/3ei(£B - £a + V), 
X + 1/3d(£b - £a - Md) - i/?ci(EB - £a + V)) Eq -J: 24) 

Ai (o,0, 0,0, 

X + i/3 p h(£B - £a) - i/3ci(£B - £a + V), 
X + i/3 P h(£B-£A)-i/3ci(£B-£A + ^)) ■ (33) 



FIG. 3: (Color Online) Renormalized phonon current from 
Eq. (15) versus bias voltage for different electron bath tem- 
peratures. At zero bias voltage (compare zoomed inset), 
bosons may either be created (higher electronic tempera- 
tures, dotted blue curve) or absorbed (lower electronic tem- 
peratures, dashed red curve). Parameters have been cho- 
sen as PuAE = /3t/(e B - £a) = 1-0, /3 D AE = 2.0 (leading 
to PphAE ~ 1.31 at symmetric tunneling rates Tu = Td), 
r L = T R = Tu = T D , /3 D U = 20, fiu = fJ,D = and 
Mi = +V/2 = -fin. 



D. Numerical Sanity Check 

We have also computed the fluctuation theorem in 
Eqns. (28), (31), and (34) by performing the required 
one- or two-dimensional integration using the dominant 
eigenvalue numerically. Within the boundaries of numer- 
ical accuracy, we have found complete agreement with 
our results (not shown). Naturally, we have also tested 
the independence on electronic tunneling rates and the 
simpler dependence on an average boson temperature. 



IV. SUMMARY 

We have investigated the multi-terminal fluctuation 
theorem for full counting statistics for a four-terminal 
model including bosonic and fermionic channels and 
Coulomb-interaction as well as phonon-assisted electron 
tunneling. We find that under the Born-, Markov-, and 
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secular approximations that under the assumption of a 
nondegenerate system spectrum lead to the conventional 
rate equations, the complete FT is fully satisfied. As 
these rate equations satisfy local detailed balance at each 
terminal separately, this result was expected. 

However, when not the complete information is gath- 
ered on all energy and matter fluxes, the FT may be ap- 
parently violated, formally expressed e.g. by a renormal- 
ized bias voltage. Thus, the modification of the FT may 
be qualitatively different from QPC back-action models 
but rather mimic the FT found for a Maxwell demon 
model. In addition, when reservoirs at different thermal 
states couple identically to the system, they may act as 



a single bath at some average thermal state - which may 
destroy universality of the FT (independence of the tun- 
neling rates). 
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